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ASYMPTOTIC LAWS FOR COMPOSITIONS DERIVED 
FROM TRANSFORMED SUBORDINATORS^ 
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Utrecht University, University of California and University of Paris VI 

A random composition of n appears when the points of a random 
closed set TZ C [0, 1] are used to separate into blocks n points sampled 
from the uniform distribution. We study the number of parts K„ of 
this composition and other related functionals under the assumption 
that TZ — <j}{S,), where {St,t > 0) is a subordinator and (p: [0, oo] — > 
[0, 1] is a diffeomorphism. We derive the asymptotics of Kn when 
the Levy measure of the subordinator is regularly varying at with 
positive index. Specializing to the case of exponential function <l>{x) = 
1 — e^^ ,we establish a connection between the asymptotics of Kn and 
the exponential functional of the subordinator. 

1. Introduction. A composition of n with positive integer parts may 
be represented by a configuration of stars separated by bars, for instance, 
x:*! X:***!* encodes the composition (3,2,4, 1) with weight 10, length 
4 and four parts 3, 2, 4, 1. A stochastic analogue of this construction appears 
when we assume the points of a closed random set TZ C [0, 1] in the role of 
bars, and n independent random points sampled from the uniform distribu- 
tion on [0, 1] in the role of stars, see [11, 13, 15, 16, 25]. Given this data, we 
define an ordered partition of the set {ui, . . . ,Un} by assigning two points 
Ui < Uj to the same block if and only if Uj and Uj are not separated by TZ, 
meaning that TZO [ui,Uj] = 0. That is to say, Uj G TZ forms a singleton block, 
and if Uj and Uj fall in the same gap (open interval component of [0, 1] \7^), 
then these points are assigned to the same block. A composition Cn is defined 
to be the record of block sizes, ordered from left to right. Exchangeability 
in the infinite sample ui,U2, - ■ ■ results in a simple consistency condition of 
Cn's as n varies, that is, the sequence (C„) is a composition structure in the 
sense of [11, 12, 16]. 
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The model just described offers a general framework for a wide range of 
"species sampling" problems, as studied in statistics and population genet- 
ics. In these applications one postulates some idealized infinite population, 
randomly partitioned into various species, with a total order on the set of 
species. A sample from such a population is understood as an exchangeable 
sequence of random variables (Xj), and a composition C.„ is defined as the 
record of multiplicities of distinct values represented among X\ , . . . , ? 
the order of increase of the values. Then (Cn) is a composition structure and 
by a de Finetti-type result [11], it can be uniquely associated with some ran- 
dom closed set TZ C [0, 1], which appears as a way to uniformize the limiting 
empirical distribution of (Xj). 

Let Kn be the length of C„ (this variable may be interpreted as the number 
of distinct species in a sample). The growth properties of moments of Kn 
are sensitive functions of the random set TZ. Logarithmic and power-like 
asymptotics of the moments are known in the case when TZ is derived by 
scaling the range of a suhordinator {St,t € [0, T]), that is, increasing process 
with stationary independent increments restricted to a finite time interval 
[1, 24, 25, 26]. (See [3] for general background on subordinators.) 

In this paper we study asymptotic properties of Kn for the random sets 
obtained by transforming the unrestricted range of a subordinator. Specif- 
ically, we consider TZ = (j){S,), where {St,t > 0) is a drift-free subordinator 
and (j): [0, oo] — > [0, 1] is a diffeomorhism. We assume the Levy measure u of 
the subordinator to be regularly varying in the sense that 

(1) u[y,^]=iil/y)y-", y [0, 

where < a < 1 and the function £ is slowly varying at oo. We also consider 
the process Kn{t), the number of parts of the partial composition produced 
by the transformed subordinator restricted to the time interval [0,t]. Other 
quantities of interest are Kn,r and Kn,r{t), defined as the multiplicity of part 
r in Cn and multiplicity in the partial composition, respectively. 

We show that, as n ^ oo, the length Kn is asymptotic to a power-like, 
regularly varying function of n multiplied by a random factor L. The factor 
L is identified explicitly as an integral functional of the subordinator. Similar 
results also hold for Kn,r, Kn{t) and Kn,r{t). The appearance of a random 
factor is due to variability in the gap sizes, as can be compared with a result 
by Karlin [20] which states that the number of distinct values in a large 
sample from arbitrary nonrandom discrete distribution is asymptotic to the 
mean number of such values. 

In the special case (piy) = 1 — e~^, the set 7^ is a multiplicative subordina- 
tor and the composition C„ inherits a characteristic regenerative property 
from this set [14, 16]. We show that L specializes in this case as the well- 
known exponential functional of subordinator. The distribution of L is then 
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uniquely determined by the power moments which are given by a known 
formula reproved here by elementary tools in the case of subordinators. 

In the regenerative case, the distribution of Kn is well known for the 
composition described by Ewens' sampling formula, in which case Kn is 
of logarithmic growth [1, 25]. More generally, Gnedin [13] has previously 
shown that the logarithmic growth of Kn is typical when the Levy measure 
is finite. For compositions belonging to the two-parameter family [16, 21, 25], 
the proper format for Kn is for parameters (a, 9) with < a < 1. Another 
interesting case is that of slow variation, when the relation (1) holds with a = 
and some i{l/y) exploding at 0. This includes the gamma subordinators 
whose Levy measure has a logarithmic singularity. This case is very different 
from the case of regular variation with positive index a and is being treated 
separately [2, 17]. 

We shall be assuming throughout that (1) holds, which entails that the 
Levy measure is infinite. When the Levy parameters {v, d) are multiplied by a 
positive factor c, the variables Kn,r remain unchanged, but Kn^r{t) should be 
replaced by Kn,r{t/c)- Basically, we assume that the Levy measure satisfies 
z^{oo} = and that the drift coefficient is 0, unless explicitly stated. 

It should be mentioned that there are many other constructions of random 
compositions, but typically these compositions are not consistent as n varies. 
One obvious possibility, in terms of the "stars and bars" representation, is 
to exploit the Bernoulli scheme, that is, to allocate a bar at each possible 
position with fixed probability p. (The particular choice p = 1/2 corresponds 
to the uniform distribution on the set of all compositions of weight n, see 
[18].) The expected length of such composition grows linearly with n, while, 
for composition structures, we have = o(n), provided the Lebesgue 
measure of 7^ is (which means that the positive frequencies of distinct 
species sum to 1). See [15] for a complete characterization of the composition 
structures obtained by truncating a single infinite sequence of stars separated 
by bars at positions visited by some increasing random process on integers. 

The rest of the paper is organized as follows. In the next three sections we 
modify Karlin's results on occupancy problems, we provide some analysis of 
the gap counts necessary to apply these results to the composition derived by 
a general transform of subordinator and we formulate the strong laws for Kn 
and the like. We specialize then to multiplicative subordinators in Section 
5. In Section 6 we continue to consider the regenerative case, but replace 
fixed-n sample by a Poisson point process, we then analyze recursions for the 
moments of the length of poissonised composition and show the convergence 
of the scaled moments of Kn- 

2. General strong laws. Karlin [20] studied the number of different types 
represented in a sample from a fixed discrete distribution with infinitely 
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many positive masses. His results open a clear path to the strong laws for 
Kn and Kn,r- Let TZ be an arbitrary closed subset of [0, 1] with zero Lebesgue 
measure. Let Cn be the composition derived from TZ by separating uniform 
points. Conditionally given TZ, the number of parts of C„ is the same as the 
number of different types represented in a sample from the discrete distri- 
bution with masses equal to the gap-sizes. Therefore, by [20], Theorem 8, as 
n — > oo, 

(2) if„~E(K„|^), Kn,r--nKn,r\n), r>l, 

where ~ means that the ratio converges to 1 almost surely. For x > 0, let Nx 
be the number of gaps of TZ of size at least x. The following is a variation 
of [20], Theorem 1, equation (23) and page 396. See also [25], Lemma 34. 

Theorem 2.1. Let i be a positive slowly varying function and L a non- 
negative random variable. The convergence 



Nx 



i{l/x)x-°' 
with < a < 1 implies, for n — oo, 



a.s.,x J, 0, 



Kn j^f^ ^r Kn,r aT{r - a) 



almost surely, and the same convergence with a = 1 implies 

^" -L, J^Ih]_^L, 4rT^ / ^ forr>l 



ni*{n) ' ni*{n) ' ni{n) r(r — 1) 

almost surely. Here i* is another function of slow variation at oo, defined 
by the converging integral 

poo g-i/y 

(3) t{t)= £ity)dy. 

Jo y 

Proof. Let us start with Kn- By (2), it is sufficient to determine the 
asymptotics of conditional expectation. To this end, we introduce a random 
measure 7 on ]0, 1] by defining its tail 

j{x):=-f[x,l]=Nx, xe]0,l], 

to be the number of gaps of TZ of size at least x. The measure 7 is atomic 
and assigns to each x € ]0, 1] an integer weight equal to the number of gaps 
of TZ of length x. For a particular gap of length x, the probability that at 
least one of n uniform sample points hits this gap is 1 — (1 — x)'", so 

(4) ¥.{Kn\TZ) = f\l-{l- xDj{dx) =n f\l- 2;)"-^7(x) dx, 

Jo Jo 
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where the second equahty is obtained by integration by parts. Observe that 
the formula 

(5) / 7(2;) dx = xj{dx) = 1 

Jo Jo 

simply says that the total length of gaps equals 1, thus, the measure 7(x) dx, 
X € [0, 1] , is a probability measure with nonincreasing density, which takes 
only nonnegative integer values. In the last integral in (4) we recognize a 
Mellin transform and standard Abel-Tauberian arguments (see Appendix) 
imply that, for < a < 1, 

7(2;) ~ x-°'i{l/x)L for X i iff 

n[ {1 - x)"'~^j{x)dx -^T{1 - a)n"i{n)L for n ^ 00 
Jo 

and the result follows in this case. In the case a = 1, the Mellin integral is 
asymptotic to the Laplace integral 

POO 

/ e-"^2;-V(l/x)dx, 
Jo 

which converges due to (5), and becomes (3) upon substituting nx = 1/y. 
The slow variation claim for £* is Lemma 4 in [20]. 
For Kn^r, we have a similar integral representation 

(6) E{Kn,r\n)=(^''^^ ^x'-il-xr-'-jidx), 

which is obtained by a formal binomial expansion of 1 — (1 — x)". The for- 
mula follows by observing that a gap of length x is hit by exactly r sample 
points with probability Q)x^{l — x)""*". A Tauberian argument applied to 
the measure x^'y{dx) yields 

E(K„,,|^)~(-lrl(^^)E(K„|7^), 

which ends the proof. □ 

Remarks. For the two slowly varying functions in the theorem, we have 
£{t)*/£(t) — > 00 as t ^ 00 (formula (13) in [20] is misprinted), that is, for 
e{t) = (logt)-",n > 1, we have £*{t) ~ (u- l)-i(logt)^-". See [6], Chapter 3, 
for results involving two slowly varying functions like £ and £* . The relation 
between asymptotics of (4) and (6) is an instance of "smooth variation" 
properties [6], Section 1.8 of the Bernstein function defined by (4). 

Under the assumption of Theorem 2.1, the conditional distribution of Kn 
given TZ approaches a normal distribution as n — >• 00, by [20], Theorem 4 
(also see [8], Theorem 2). Karlin's results also imply a multivariate normal 
limit for the conditional distribution of the sequence {Kn,r,f ^ 1)- 
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3. Counting the gaps. Let {St,t > 0) be a subordinator with the drift 
coefficient d = and a Levy measure u satisfying i/{oo} = 0. The jumps of 
(St) correspond to the gaps of 7^ C [0,oo], which is a topological Cantor set 
provided the Levy measure is infinite. 

Let (f): [0, oo[— > [0, 1[ be a diffeomorphism, that is, a continuously differ- 
entiable function satisfying i;^(0) = 0, (/>(oo— ) = 1 and (p'it) > 0. For TZ := 
0(7^) C [0, 1], the gaps comprising TZ^ = [0, 1] \ 7^ correspond to the jumps 
of the subordinator transformed by (p. Let Nx{t) be the number of jumps 
of size at least x for the transformed subordinator restricted to [0, t] , and 
let = Nx{oo) be the number of such gaps of TZ without restriction. We 
are interested in the asymptotics of these gap counts for small x. A similar 
analysis has appeared in [23] in the case of stable subordinators. 

The analogous question for the original subordinator (St) is easy. Let 
Ny{t) be the number of gaps of TZ of size at least y, generated by the subor- 
dinator restricted to [0,t]. The counting process {Ny{t),t > 0) is a Poisson 
process with rate I'iy), thus, for small y, the behavior of this process is ruled 
by the strong law of large numbers: 

(7) Ny{t) u{y)t, yiO, 

almost surely for all t. We shall see that translating this behavior into similar 
results for Nx{t) and Nx = Nx{oo) amounts to a change of variable formula 
which was stated in [23], albeit under different assumptions on (p. 

Speaking more broadly, we may wonder about the conditions on (j) and z/ 
which imply an asymptotic relation analogous to (7) of the type 

(8) Nxit)r^i;ix)L{t), xlO, 

where tp is a scaling function and {L{t),t G [0,00]) is a positive random 
process. A principal new effect appearing in (8), as compared to (7), is that 
a nonlinear transformation of subordinator leads to a genuinely random 
scaling limit. The next question to ask is whether such a relation holds with 
some L for t = 00 and whether L(oo— ) = L, and we shall find conditions 
when this is true. 

For a Levy measure as in (1), we shall use the scaling function 

ip{x) = x-°'e{i/x). 

3.1. Finite t formula. Let Ny{ti,t2) be the number of jumps of {St,t € 
[^15^2]) of size at least y, with the convention that this is zero for ti > t2- 

Theorem 3.1. // the Levy measure satisfies (1), then for each diffeo- 
morphism 4>: [0, oo[^ [0, 1[ and <t<oo, the convergence, as x — > 0, 

(9) Nx{t)/i;{x)^ f\4>'{Su)rdu 

Jo 
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holds in the mean for each t, as well as almost surely, uniformly in t hounded 
away from oo. 

Proof. Consider a partition of [0, t] by points = to < < • ■ • < ^fc = * j 
with probability 1, each tj is a continuity point of the subordinator. As is 
easily seen, 

^x/4>'{s,,^)i'tj,tj+i) < Nx{tj,tj+i) < iV^./<^,(5^^)(tj,tj+i), 

where and rjj are the points where (j)' attains the maximum and the 
minimum on [Stj,Stj^-^], respectively. Taken together with (7), this implies 

(10) {t,+i-t,)i;(^^^^ «iv,(t^.,t.^,)«(t^.^,_t^.)^(^_|_^, 

where the notation X <^Y for positive random quantities depending on x 
means that esssupX/y < 1 for x | 0. From this and the assumption on v, 

k—l k—1 

V^(x) ^(<^'(5,^,))"(t,+i - t,) « N,{t) « iP{x)Y,{<p'{S^^)nt,+, - t,). 

j=0 j=0 

We see that Nx(t), for x | 0, can be squeezed between an upper and a lower 
Riemann sum; thus, sending the diameter of partition to zero and using the 
continuity, we obtain the almost sure convergence (9). Using the obvious 
bound Nx{t) < N^j ^^y^^i{t) where the maximum is taken over [0,t], the 
convergence in mean follows by dominated convergence. □ 

There is a minor generalization of the formula for integrals with random 
upper bound. By stopping time r, we understand a random variable tak- 
ing values in [0, oo] and such that ((r A t),t > 0) is adapted to the natural 
filtration of the subordinator. 

Corollary 3.2. If the Levy measure satisfies (1), then for each diffeo- 
morphism (j): [0,oo[— > [0, 1[ with sup(p' < oo, and each stopping time r with 
Er < oo, the convergence, as x — > 0, 

iv,(T)/^(x)^ r{<p'{Su)rdu 

Jo 

holds almost surely and in the mean. 

Proof. The almost sure convergence follows from Theorem 3.1. The 
convergence in the mean is a consequence of 

lim limsupEA'a;(t, r)/'0(x) =0, 
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which, in turn, follows from 



( 



sup <p' 



X 




by application of Wald's identity. 



ENy{T) = u{y)ET, 



and the fact that E(r — — > for t — > oo. □ 

The corollary can be applied to a subordinator killed at an independent 
time r. For example, when z^{oo} > 0, the subordinator jumps to infinity at 
an independent exponential time. 

3.2. Full range formula. We turn next to finding some conditions for the 
convergence 



JO 

in which case (8) holds for t = oo with L = limt^oo L{t) given by the integral 
in (11). One condition which seems very natural is the integr ability 



Granted this integrability condition only, we failed to prove or disprove 
whether the convergence holds in full generality, nor is there an obvious 
sufficient condition which would cover the cases of interest including slowly 
varying functions | logy| or 1/| logy| and diffeomorphisms with exponentially 
decaying or power-like tails. 

To secure the convergence, we shall make some additional assumptions 
about (j) and ly. The analysis is largely simplified by further assuming that 
the derivative (j)' is a decreasing function; in this case, in (10), we can set 

= tj and r]j = tj+i, and then, for any ti <t2, conditioning on St^ yields 
the inequality 



which is valid for each x > 0. 

Theorem 3.3. Assume that the Levy measure is as in (1), the diffeo- 
morphism (j) has decreasing derivative, the integrability condition (12) holds, 
and one of the following single or composite conditions is satisfied: 



(11) 




(12) 




(13) 




(i) for some constants a > 0, C > 0, the inequality £{u) < Ci{v) holds 
for u < av, provided u and v are sufficiently large; 
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(ii) there exist functions q and r such that: 

(iia) q{y) = oiiljiy)) as y [ 0; 

(iib) for some constants a > 0, C > 0, the inequalities r(l/v)v <u < 
av imply (.{u) < Ci{v) for all sufficiently large u, v; 

(iic) i'{x/r{x))(t)^{l-q{y))=o{i'{y)) as y [{); 

(iii) £ is hounded away from and oo on every compact subset of [0,C!o[ 
and a stronger integrability condition holds, with a in (12) replaced by a — 6, 
for some 5 > 0; 

(iv) the same integrability condition holds as in (iii) and there exists a 
function q which satisfies (iia), as well as — q{y)) = o{ip{y)) for y J, 0. 

Then the convergence (11) holds almost surely and in the mean. 

Proof. In view of liminf Nx/ip{x) > L{oo—), it is sufficient to establish 
the convergence of expectations 

~ /-oo 

hm ENx/ip{x) = E / {(p'iSt))" dt. 
Thanks to both convergence results, 



rt roc 

limE / {(l)'{Su)Tdu = K / {(l)'{Su)Tdu 
Jo Jo 



< oo, 



limEiV^(i)/V'(x) =E f\(f>'{Su)Tdu<oo, 

x^O Jq 

and in view of Nx = Nx{0, t) + Nx{t, oo), we only need to show that 
lim limsupEA'a;(t, oo)/V'(2;) =0. 

t-^oo x-^O 

Using monotonicity, 

oo oo y \ 

ENx{t, oo) = Y^ ENxit + j,t + j + l)<Y^ Eu( I ) 
U J^o \^(St+j)J 

oo 

= EY,{4>'{St+j)ri{<P'{St+j)/x)x-'-, 

j=0 

whence 

X °£(i/x) \ e{i/x) 

The rest of the proof amounts to estimating the right-hand side of this 
formula. 
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Let Tx be the first passage time over the level (<^')^(x), when (j)'{St) drops 
below X. As a consequence of Corollary 3.2, the sum of terms in (14) with 
t + j<Ta is negligible, for each fixed a > 0. And assuming (i), we estimate 
the contribution of terms in (14) with t + j > Ta by 



CE / {4>'{Su)rdu, 
Jt-i 

which vanishes for t — cxd due to (12). 

Assume (ii). Let cTz be the first passage time over the level 0^(z). By 
(iia), the contribution of terms with t + j> (Ji-q^x) is negligible because at 
most q{x)/x gaps longer than x can fit into an interval of size q{x). The 
contribution of terms with t + j < Tj.(^x) ^^^o vanishes, as x | and then 
i ^ oo, for the same reason as under the condition (i). Finally, 

ENx{Tx,ax) < e(^z7(^^^^^^ {ax - r^)+) 

by virtue of (iic), because the expected time for the subordinator to pass a 
high level y = 0^(1 — q{x)) can be estimated from the above by by, with 
some positive constant b. 

Suppose (iii) holds. By Potter's Theorem A.3(ii), there exists C > 1 such 
that £{u)/i{v) < C{u/v)~^ for all u<v. We apply this with u = 4)'{St)/x 
and V = \/x and then use the integrability with the exponent a — 5. 

Under assumption (iv) , we make use of another part of the same Theorem 
A.3(i), which guarantees the same inequality for sufficiently large parame- 
ters, say, u> A, v > A. For (St+j)) with t + j< tax we have then the 
same inequality as in (iii), thus, the contribution of such terms can be esti- 
mated as in the case (i), but with the exponent a — b. The remaining sum 
is bounded from above by 

v[A)¥.aY-q(x) 
analogously to the case (ii). □ 

The integrability condition (12) can be ensured by means of the following 
lemma found in [3], page 28. 

Lemma 3.4. For each decreasing "positive function g on [0, oo[ and each 
subordinator {St,t > 0), the following properties are equivalent: 

poo POO 

/ g{t)dt <oo <^=^ / g{St)dt <oo a.s. 
Jo Jo 

POO 

E / g{St) dt < oo. 
Jo 
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Remarks. The function i{l/y) = jlogy|'',/9 > 0, is decreasing, thus, 
part (i) of the theorem apphes. For cl){y) = 1 — (y + the integrabil- 

ity condition is fulfilled if a(/3 + 1) > 1, thus, selecting b in the range 1 — a < 
b < a/3, part (iv) applies for any £, with q{y) = y^. Part (ii) is useful for 
i{l/y) = I logyl^'^, /9 > 0, in which case we can take r{y) = y'' ,0 < b < 1, to 
meet (iib). 

Note that, for q{y) = y^ ,0 < b < 1, the condition on (j)^ can be refor- 
mulated in terms of (p, for example, — x) = o(x~") is equivalent to 
l-0(y)=o(y-i/"). 

4. Strong laws. Strong laws for Kn and the like for composition derived 
from a transformed subordinator follow by combining the results in the two 
previous sections. Introduce 

(15) L{t)= f\^'{Su)rdu, L = L(oo). 

Jo 

Recall that notation Kn{t) or Kn,r{t) refers to the parts of the partial com- 
position produced by the range of {(j){Su),u G [0,t]). 

Theorem 4.1. For < a < 1, the regular variation assumption (1) im- 
plies 

Knit) Kn,r{t) 



r(l - a)n'^e{n) 



1(1 — ajn"i;(nj V / 



with probability 1, as oo. And if (p satisfies also the conditions of The- 
orem 3.3, we have 



(-ir^(") 



r(l-a)n"^(n) ' r(l - a)n"^(n 
For a = 1, the analogous results are read from Theorem 2.1. 

Asymptotics in the drift case. We sketch the extension to the case of 
subordinator with positive drift d. In this case the length of composition 
satisfies Kn ~ Kn,i ~ nX(JZ), where A denotes the Lebesgue measure, thus, 
the asymptotics follow from the next lemma which generalizes [16], Corollary 
5.8. 

Lemma 4.2. Let (St) be a subordinator with drift d > 0, and cj): [0, oo] — > 
[0, 1] be an absolutely continuous strictly increasing function. Then 

A(0(7^)) = d/ (t>'iSt)dt. 
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Proof. Because (St) is almost everywhere differentiable with derivative 
d, by the change of variable St = y, 

d ct>'{St)dt= / <i>'{y)dy = \{m))- □ 

In the drift case it is sensible to distinguish between genuine single- 
ton parts which are caused when a sample point hits 7^, and the other 
occasional singletons induced by open gaps comprising Tf^. (For fixed n, 
the composition C.„ does not allow one to make this distinction.) Denoting 
Kn^i- and i^n,i+ the counts of genuine and occasional singletons, we have 
Kn = Kn,i- + Kn,i+ + J2'?^2 ^n,r, and Kn ~ -f^n,!-- The asymptotic behav- 
ior of variables i^n,i+ and i^n.r for r > 1 is then still as in Theorem 2.1, as 
follows by noting that the gap-counting Theorems 3.1 and 3.3 and Corollary 
3.2 are also valid for subordinators with drift. 

A curious phenomenon occurs for a = 1: normalizing the Levy data (z/, d) 
so that d = 1, we have all three variables Kn/n, Kn^i-/n and Kn,i+ / {ni* {n)) 
approaching the same limit. 



5. Regenerative compositions. We shall specialize the results of the pre- 
vious section to the regenerative compositions appearing when TZ is the range 
of the multiplicative subordinator 5^ = 1 — exp(—St)- In this case there is a 
simple connection between 5^ and L{t) and a nice formula for the moments 
of L. 

With each multiplicative subordinator, we associate the area process 

rt 



A{t)= I {l-Su)du 
Jo 



and its terminal value A = ^(oo) obtained by taking the infinite integration 
bound. In terms of the (additive) subordinator. 



: / ex.p{—Su)du= / {1 — Su)du 
Jo Jo 



is the widely studied exponential functional, see [4, 5, 7, 28, 29]. 

Let be the measure on [0, 1] obtained by the exponential transform 
(p{y) = 1 — e~y from the measure ly. The Laplace exponent is thus given by 

/■oo rl 

$(s)=/ {I - e-'y)v{dy) = {l-{l-xY)u{dx). 
Jo Jo 

Because (/>'(0) = 1, the assumption (1) implies that the tail of the trans- 
formed measure satisfies v[x, 1] ~ x~"l{l/x) for x | 0. For arbitrary a > 0, 
the process 5^"^ := 1 — (1 — St)°' is itself a multiplicative subordinator with 
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Levy measure i?*-"^ related to v by v^'^^x, 1] = t'[l — (1 — xY^°' , 1]. The rela- 
tion between the corresponding Laplace exponents is <^(°)(s) = <^(as). That 
is to say, 

s'f'^ = 1 — exp(— a5f ) 
is the multiplicative counterpart of the scaled subordinator {aSt). The area 
process for {s[°''^) is 

L(°)(t) := [\l-S^^^)du= l\l-SuTdu= f\xp{-aSu)du 
Jo Jo Jo 

and we define L(") =L(°)(oo) to be the A-functional for {S'^"^). 
For the scaling function iIj{x) = x~°'£{l/x), we have the following: 

Theorem 5.1. Suppose the Levy measure fulfills (1), then, for x [0, 
the jump counts of the multiplicative subordinator St = l — exp( — St) satisfy 

iV,(t)/V^(x) ^ L(")(t), N^mx) ^ 
for X 10 almost surely and in the mean. 

Proof. We have (j)'{y) = e~^. Then part (iv) of Theorem 3.3 applies, 
because the required integrability holds for any arbitrary positive power 
and the second condition is fulfilled with q(x) = x. □ 

The application of Theorem 2.1 results in the following: 

Corollary 5.2. // the Levy measure fulfills (1), then, for < a < 1, 
the composition induced by the multiplicatively regenerative set IZ satisfies, 
for n oo, 

Jf„/(n"£H) ^ r(l - a)L(°) 
almost surely and in the mean. And for a = 1, 

Kn/{nt{n))^L^^'>=A 
almost surely and in the mean, with i* as in (3). 

Generalizations to Kn{t) and Kn,r{t) follow in the same way. 

The distribution of L^") admits some exponential moments, hence, it is 
uniquely determined by its integer moments. They are given by the following 
formula which was recorded for general Levy processes in [7], though it can 
be traced back in special cases to much earlier literature (see, e.g., [10], page 
283): 



(16) 



1,2, 
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Example. Consider the two-parameter family of regenerative composi- 
tion structures, as in [16], with 

i>[x,l]=x-°'{l-xf, 0<a<l,e>0. 

In this case we have 

sTil-ans + ei 
^ ' T{s + l-a + 9) ' 

thus, 

Tie + l)(a + 9){2a + B)--- {{k - l)a + 0) 



r(A:a + 0)Q* 



in agreement with formula (192) from [25]. This specializes for 9 = to the 
integer moments of the Mittag-Leffler distribution with parameter a. 

Now we shall give a new proof of (16) in the case of subordinators. The 
method we use here is not applicable to the exponential functionals of more 
general Levy processes, as considered in [5, 7], but it is much more elemen- 
tary and apparently more natural in the context of multiplicatively regen- 
erative sets. By the above discussion, it is sufficient to prove the formula 
for a = 1, that is, for the area functional of a multiplicative subordinator. 
Letting = KA^, we wish to show that 

A;' 

(17) mfc= . " , A; = 0,1,.... 

Finite Levy measure, no drift. Suppose first that d = and z? is a proba- 
bility measure. Let Xj be a sample from u. Then (St) is a step function whose 
range is a stick-breaking sequence Xj{l — Xi) • • ■ (1 — Xj-i),j = 1, 2, ... , 
complemented by and 1. The jumps of (St) occur at the epochs of an 
independent homogeneous Poisson process. Therefore, A is representable as 
a random series 

A = EiXi + {Ei+E2)X2{l-Xi) + (^1 +^2 + ^3)^3(1 -^i)(l -^2) + • • • , 

where Ej are jointly independent exponential random variables with mean 
1, also independent from the Xj^s. The series is finite only if ?{1} > 0. We 
can also re-arrange terms and write A in the form 

A = Ei+E2{l-Xi) + E^{l-Xi){l-X2) + ---, 

from which we deduce 

(18) A = E + {l-X)A', 
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where A' is a replica of A, independent of {E,X) = {Ei,Xi). By virtue of 
the formula 

E{l-X)'' = l-^{k), 

the expectation is computed as MA = 1/$(1), in accord with the k = l case 
of (17). Furthermore, the identity 

k 

j=Q ' 

implies that the moments (rrijt) of A exist and satisfy the recursion 

mk^{k) = ^ — (1 - ^{j))mj. 
i=o ^■ 

To solve the recursion, split out the last term and substitute the same iden- 
tity but with /c — 1 to arrive at 

mk^{k) = k{l - <t{k - l))mfc_i + kmk^i^{k - 1), 

which is the same as the simple multiplicative recursion 

mk^{k) = knik^i, 

whose unique solution with the initial value niQ = 1 is that given by (17). 

Replacing the probability measure ? by its positive multiple, say, i^c = cu, 
we can write a series representation for the corresponding functional Ac 
exactly as above, but with Ej/c instead of Ej. Since E(i?/c)^ = k\/c^, the 
same computation yields the additional factor in the denominator. This 
agrees with (17) because the new Laplace exponent is the multiple c^. 

Finite Levy measure, positive drift. The moments formula for a subordi- 
nator with drift can be proved analytically using approximation by drift-free 
subordinators, but it is more instructive to inquire into this case separately. 
Assuming u[0, 1] = 1 and d > 0, the subordinator coincides with function 
1 — e"'^* for t < El, and the jump at time Ei is e~'^^^Xi. The first-jump 
decomposition is 

A £ e-^^{l - X)A' + (1 - e-^^)/d, 

with the same convention as in (18). The recursion for moments is obtained 
by using 



E{1 - X)'' = 1 - <^>{k) + kd 
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and exploiting the fact that e is distributed according to Beta(l/d, 1), 
we obtain 



The solution is again (17), as justified by the same inductive argument. 
Repeating the above scaling argument, we see that the formula also holds 
in the case d > and arbitrary finite v. 

General suhordinator. Given arbitrary Levy data (d, consider the 
family of subordinators with parameters (d, z/^), where is a truncated 
measure that coincides with v outside [0,e] and is zero within [0,e]. Using a 
version of the well-known recipe, we can construct the corresponding multi- 
plicative subordinators St and Se^t using the same Poisson point process in 
the strip [0, cxd[ x [0, 1] with intensity measure Lebesgue x z?, so that 



where the product is over the atoms (rj,Aj), and S^^t has a similar repre- 
sentation, with the only distinction that the factors corresponding to the 
atoms with Aj < e do not enter into the product. By construction, S^^t T St 
as e I and the convergence is uniform in t. Thus, by monotone conver- 
gence, we have for the corresponding integrals yle(ti,t2) i ^(^11^2) almost 
surely and with all moments, for all < ti < t2 ^ 00. But all measures Vgr 
are finite, thus, as we have shown, the moments formula is true, therefore, 
the formula also holds for (d, u) since the corresponding Laplace exponents 
satisfy T 

Remark. For A' a copy of A independent of {St-,A{t)) (fixed t), we have 



^(fc) 
kd + l 




fc! (l-«^(j)+jd) r(i + l/d) 
j! d'^-i+i r(yt+l + l/d) 



A{t,oo)^{l-St)A', 



A ^ A{t) + (1 - St)A'. 



This leads to 



¥.{A{t,oo)f 



exp(-t$(A;)) 



¥.A{t) 



1 -exp(-t$(l)) 
W) 



Higher moments of A{t) are not immediate because of dependence between 
St and A{t). See [7] for the formulas with t replaced by an independent 
exponential variable. 
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6. Poissonized compositions. A closely related type of structure appears 
when the uniform sample of fixed size n is replaced by a Poisson point pro- 
cess of rate p. A composition of random weight n{p) appears by separating 
the Poisson points into blocks by means of a random closed set TZ C [0, 1]. 
We shall denote this poissonized composition Cp, and provide with all 
quantities related to it. The relation to the fixed-n composition is therefore 
Cp = Cn, conditionally, given n{p) =n. 

Poissonization is useful for two reasons. Generally speaking, it is a pow- 
erful technique for asymptotic considerations in combinatorial problems, 
allowing one to explicitly exploit the spatial independence where other- 
wise only a kind of asymptotic independence is available (see, e.g., [27] for 
overview). On the other hand, poissonization yields a family of compositions 
{Cp,p > 0) which satisfies a consistency condition analogous to the defining 
property of partition or composition structures [11, 16, 25]. Explicitly, for 
any p> and x G [0,1], probability distributions of the following compo- 
sitions coincide: (i) a composition with rate parameter p{l — x) and (ii) a 
thinned composition which appears when the atoms making up a sample 
with rate p are deleted independently with probability x. 

In the sequel we shall only consider the case when TZ is the range of a 
multiplicative subordinator St = 1 — exp(— 5^), in which case the distribution 
of (i) also coincides with the distribution of (iii), a tail composition of the 
composition of rate p which appears to the right of x, conditionally, given 
X € TZ. The last equivalence is analogous to the regenerative property of 
the fixed-n compositions [16]. The same composition of random integer n{p) 
appears when the range of additive subordinator (St) is used to separate into 
blocks atoms of inhomogeneous Poisson process on [0, oo] with exponential 
intensity measure pe~^ dx,x G [0,oo]. We denote by Kp the length of the 
poissonized composition, and Kp{t), Kp^r{t) stand for the number of parts of 
the partial composition produced by the range of multiplicative subordinator 
up to time t. 

We proceed with the convergence results which recover and complement 
the results in the previous sections. The equivalence of strong laws for 
(Cp) and (C„) is quite obvious, and for quantities like moments, there is a 
well-developed analytical technique of poissonization/depoissonization [19], 
though we shall use more elementary arguments. 

6.1. Recursions. Let J-t be the a-algebra generated by the subordinator 
{Su,u G [0,^]) and by the Poisson configuration on [0,54]. By the indepen- 
dence property of the Poisson process^ and the regenerative property of TZ, 
the tail composition induced hy TZri]St,l] is independent of Tt. This obser- 
vation is a source of recursions related to the poissonized composition. 
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Let Pj{p),j = 0, 1, . . . , be the distribution of Kp. Each pj may be extended 
to an entire function of the complex variable, with initial value Pj{0) = 
[with the only exception po{p) = e~^]. Introduce the factorial moments 
f^){p)=tKp{Kp - 1) . . . (Kp - m + 1), with m = 0, 1, . . . , /^(p) = 1. 

Lemma 6.1. The following integral recursions hold. For j = 1, 2, . . . , 

^\p,{p)-e-P-p,{p{'^-x))D{dx)) 



(19) 



\l-e-nVj-i{p{l-x))V{dx), 



and for ?7t, = 1, 2, . . . , 

\f^'^\p)-f^'^\p{l-x)))v{dx) 



(20) 

= m 



[\l-e-nf^"'~'Hp{^-x)Mdx). 
Jo 



Proof. Each Pj{p) may be written as a generating function whose co- 
efficients are rational functions in the variables ^{n),n > 0, for example, the 
probability of one-part composition is 

„ ~p"$(n:n) 

Pi{p) = e 



n=l ^ ' 



where 



$(n:m)= (^^^ ^\™(1 - x)"J?(dx), 1 < m < n. 

Thus, the statement is of purely algebraic nature and can be translated as 
a series of polynomial identities in these variables. Thus, it is sufficient to 
consider the "stick-breaking case" of finite Levy measure, normalized to a 
probability measure, when the recursion is proved by conditioning on the 
first break X = x with distribution z7. Indeed, there are j blocks when either 
[0,x] contains a Poisson atom and then [x,l] generates j — 1 blocks [with 
probability pj^i{p{l — x))], or [0,x] is empty and [x,l] generates j blocks. 
This gives 

pM = f\e~P-p,{pa - x)) + (1 - e-^Pj-iiPil - x))Mdx). 
Jo 

To keep this formula right for arbitrary finite v, we should put pj (p) into the 
integral, then the formula becomes homogeneous in the $(n)'s and holds in 
general, for algebraic reasons. 
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To prove (20), start with the definition 

oo 

/M(p)=^p^.(p)j(j-l)...(j-m + l), 

i=o 

then multiply both parts in (19) by j{j — 1) • • • (j — m + 1) and sum over j 
using the identity 

j(i - 1) • • • (j - "1 + 1) = (j - 1) ■ • ■ (i - m) + m{j - 1) ■ • ■ (j - m + 1). □ 

Manipulation with power series allows, in principle, computing the distri- 
bution of Kn with all moments, for fixed-n compositions. Let us demonstrate 
this on the expectation f^^\p) = ¥.Kp. For fn^ = EX^, we have the pois- 
sonization identity 

n=0 

Substituting into (20) and integrating, we obtain a relation between gener- 
ating functions 

OIL. I i. 

n=l 

Multiplying by and extracting the coefficients, we get 

which is a familiar expression for KKn, see [16]. 

6.2. Asymptotics. The convergence ^p/(r(l - a)p''e{p)) L(") a.s. fol- 
lows exactly as in Sections 4 and 5 for < a < 1 (and with a proper scaling, 
also for a = 1), in the footprints of Karlin [20], where the Poisson model was 
treated in parallel with the fixed-n case. In this section we show that the 
recursion (20) implies the moments formulae analogous to (16). This can be 
regarded as a proof that the convergence holds together with all moments, 
and also as yet another derivation of the moments formula (16). 

Introducing the poissonized Laplace exponent (not to be confused with $ 
written in terms of u) 

$ip)= [\l-e-nHdx), 
Jo 

we have as p ^ cx) 

g ( r{l- a) p^£{p), forO<a<l, 
^^'^\pr{p), fora = l, 

(see the Appendix). 
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Proposition 6.2. The factorial moments satisfy 

with c^"^^ given by 

for < a < 1. For a = 1, the factors T{1 — a) should be omitted and i 
replaced by £* . 

Proof. We concentrate on the case < a < 1, leaving the case a = 1 to 
the reader. Trivially, f^'^\p) = 1- Suppose by induction that the asymptotics 
holds for some m — 1. Then setting b = T{1 — a) and g{p) = {p°'i{p))'^ , we 
have 

(22) [\l - e-^^)/(™-i)(/>(l - x))V{dx) ~ 6c("^-^)5(/3). 

Jo 

This is shown by splitting the integral at e and replacing the integrand for 
X G]0,e] by its asymptotics. To justify the induction step, fix e and suppose 
there exists arbitrarily large p such that 

/(-)(p)>(l + e)c(™)5(p) 

(we wish to lead this assumption to a contradiction). Then, perhaps selecting 
e smaller, for any fixed constant C, there exists arbitrarily large p such that 

/(-)(p)>(l+e)c(-)5(p) + C 

[just because f^^\p) oo]. Up to the end of this paragraph, p = p{C) will 
be the minimal p for which the inequality holds. Note that p{C) — > oo as 
C — > oo. Thus, we have 

/(™)(P) = (1+6)C(™)5(P)+C, 

/("^)(p(l - x)) < (1 + e)c('")5(p)(l - x)'^" + C, X G ]0, 1], 

and substituting this into (20), we see that the left-hand side is estimated 
from below by 

(1 + e)c^"'^g{p) [\l - (1 - x)''"')D{dx) = (1 + e)c('")g(p)$(am) 
Jo 

= (l + e)c(— i)65(p), 

where /? — > oo and we used monotonicity. This disagrees with the right-hand 
side of (20) given by (22), giving the required contradiction. Thus, 

/M(p) 
lim sup . . , , < 1 + e. 
c(-)(7(/9) 
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A symmetric argument shows that 

liminf , , >l-£. 

Letting e — > ends the proof. □ 

Depoissonization follows rather easily. Recall that the collection of atoms 
of the Poisson process with rate p can be identified with a uniform sam- 
ple {lii, . . . , n„(p)}, with n{p) distributed according to Poisson [p). By the 
obvious monotonicity of we have 

Knl{n{p{l - e)) <n)> Kp(i„,)l(n(p(l - e)) < n), 

Knl{n{p{l + e)) >n)< K p^i+,)l{n{p{l + e)) > n). 

Selecting n = p and letting p^ oo, the elementary large deviation bounds 
for the probability P(/9(l — < n{p) < p{l + e)) imply that, for n oo, 

Kn{l~e) <.Kn<. 

Letting e | and using Proposition 6.2, we see that Kn ~ Kn almost surely 
and with all moments. 

Observing that the computation of the constants (21) is equivalent to 
the formula (16) [thus, we have yet another proof for (16)], and recalling 
Theorem 2.1, we summarize the above discussion in the following theorem. 

Theorem 6.3. The almost-sure convergence 

^p/(p"^(p))-r(l-a)L(-), p^oo, 

Kn/in^lin)) ^ r(l - a)L(") , n ^ oo, 

holds for both Poisson andfixed-n compositions together with the convergence 
of all integer moments for < a < 1, and with a proper scaling also for a = l. 

The convergence of moments of Kn,r, Kn{t), Kn,r{t), Kp^r, Kp{t), Kp^r{t) 
(with obvious definition of the last two random variables) follows from the 
theorem by dominated convergence. 

6.3. A martingale approach. Extending the discussion in the previous 
section, consider Kp{t), which is the number of blocks of a poissonized com- 
position produced by the subordinator up to time t. We can view {Kp{t),t > 
0) as either an increasing process with unit jumps or a point process of those 
jump-times of (St) which have jump intervals covering some sample points. 
The compensator for Kp{t) is 

Ct= f^p{l-Su))du. 
Jo 
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By observing that 1 — x)) is the probabihty that a gap with leftpoint 
X is hit by a Poisson atom, the formula can be first argued in the renewal 
case. The general case follows by extrapolation from the case of finite Levy 
measure. This readily implies the following: 

Lemma 6.4. For each p> 0, the process 

Mt:=Kp{t)-Ct, tG[0,oo], 

is a square-integrable martingale with unit jumps and quadratic predictable 
characteristics {M)t = Ct- Furthermore, 

EM^ = E(Kp{t) - Ctf = ECf. 

Proof. The squared jump magnitudes of Mf are 1. This implies that the 
submartingale Mf has the same compensator as Kp{t), that is, {M)t = Ct, 
as in [22], Section 6.2. □ 

The lemma opens yet another approach to the convergence results, for 
which we give below the jC^-version. Note that the scaling by $(p) is asymp- 
totically the same as the scaling by ^{p) (see Appendix), which is asymptotic 
to r(l - a)i{p)p'' for < a < 1 and to e*{p)p for a = 1. 

Theorem 6.5. Under the regular variation assumption (1), as p—>oo, 

^ ^ Hp) 

almost surely and in C? . An analogous result is valid for Kp{t) for each 
t>0. 

Proof. We wish to establish the convergence (23) in C^. Use Lemma 
6.4 to obtain 

(24) e(3^- rm^,ty=j-Erm^,t. 

V$(p) Jo $(p) / $(p) Jo <5(p) 

Also observe that 

for,-, $(p(l-gf)) ~ 

(25) ^— >{l-St) asp^oo 

$(p) 

almost surely for each fixed t. Thus, (23) would follow by dominated con- 
vergence once we could bound 
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from above by a square- integrable random variable. To this end, write ^{p) = 
p"-^o(p) with slowly varying Iq [so £o{p) ^ i{p)], then, if the Potter's bound 
were valid for £q on ]0, cxd[, we could estimate 

^p{i-St))Mp)<c{i-Str-^ 

with some small 6 > and C > 0, whence 

roo ^ roc ^ 

/ (<I>(p(l - St))/Hp)) dt<C (1 - StT-^ dt e C\ 
Jo Jo 

To make this argument precise, we fix some sufficiently large constant c = X{C, 6) 
required in Theorem A.3(i), and then split the integral at the first passage 
time (Ji_c/p of (St) over the level 1 — c/p. The tail integral 




^p{l-St))dt 



is bounded from above by a Poisson random variable with mean c, which is 
obviously square-integrable, for each c > 0. Thus, it is sufficient to exploit 
Potter's bound for Iq on [c, oo[. □ 



APPENDIX 

A.l. Abel-Tauberian theorems. An exposition of Abel-Tauberian the- 
orems for the Laplace transform is given in [9], Section XIII. 5. Bingham, 
Goldie and Teugels ([6], Section 1.7 and Chapter 4) give a fuller account, 
also for more general integral transforms, including the Mellin transform. 

We establish next some elementary connections between the integral trans- 
forms in a form suitable for applications to subordinators. Consider the two 
Laplace exponents (also called Bernstein functions) 

poo pi 

$(s):=/ {I - e~'y)v{dy) = {I - {I - xy)u{dx) 
Jo Jo 

and 

p\ poo _ 

$(s):=/ {I - e-''')u{dx) = {l-e-'^^-^~'^)v{dy), 
Jo Jo 

where the measure i'[dx) on ]0, 1] is the image of i^(dy) on ]0, oo] via x = 
1 — e~^, and it is assumed that ^>(s) < oo for some (and, hence, all) s > 0. 
The function $ is the poissonization of $, that is. 



oo n 

1 1. 
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Lemma A.l. Whatever the Laplace exponent ^, 

lim (^>(s) -$(s)) = 0. 

s — ^oo 

Proof. With a hint from [8], page 1257, we have, for < x < 1 and 
s>0, 

< e-'"" - (1 - xY < sx^e-"^ < e'^x. 

Using this we get 

< $(s) - = /\e-*^ - (1 - xy)u{dx) < e-^ xu{dx) = e"^$(l), 
Jo ^0 

so the claim follows by dominated convergence. □ 

Corollary A. 2. Ahel-Tauherian relations (as s ^ oo) between the tail 
of the measure i/[l/s,oo] and the Laplace exponent are the same for $(s) 
and 

Estimates of the difference can be given under the assumption of regular 
variation. The Laplacian case of monotone density, relating asymptotics of 

and <1>, is covered by a combination of Theorems 3 and 4 in [9], Section 
13.5. 

A.2. Potter's theorem. 

Theorem A. 3 ([6], Theorem 1.5.6). Leti be a function of slow variation 
at infinity. 

(i) For arbitrarily chosen constants A> 1,6 > 0, there exists X = X(A, 6) 
such that 

£{y)/£{x) < Ama^{{y/xY, {y/xY^) {x>X,y> X). 

(ii) Lfi is bounded away from and oo on every compact subset of [0, cx3[, 
then, for every 6 > 0, there exists A' = A' (6) > 1 such that 

i{y)/i{x) <A'max{{y/xY,{y/x)-^) (x>0,y>0). 
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